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Abstract 

Partial eigensystem assignment with output feedback can lead to an unstable closed-loop 
system. However, output feedback with passive linear time-invariant systems, such as flexible 
space structures, is guaranteed to be stable if the controller is dissipative. This paper presents a 
novel approach for synthesis of dissipative output feedback gain matrices for assigning a selected 
number of closed-loop poles. Dissipativity of a gain matrix is known to be equivalent to positive 
semidefiniteness of the symmetric part of the matrix. A sequential procedure is presented to 
assign one self-conjugate pair of closed-loop eigenvalues at each step using dissipative output 
feedback gain matrices, while ensuring that the eigenvalues assigned in the previous steps are 
not disturbed. The problem of assigning one closed-loop pair is reduced to a constrained solution 
of a system of quadratic equations, and necessary and sufficient conditions for the existence of 
a solution are presented. A minimax approach is presented for determining parameters which 
satisfy these conditions. This method can assign as many closed-loop system poles as the 
number of control inputs. A numerical example of damping enhancement for a flexible structure 
is presented to demonstrate the approach. 
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Introduction 


Pole placement has been studied extensively in the literature to shape transient response 
of the closed-system. Many techniques are available for eigensystem assignment with full 
state feedback. These techniques exploit the freedom beyond pole placement for multivariable 
systems, identified in Ref. 1, for other objectives such as enhancing robustness of the pole 
placement and minimizing the norm of the feedback gain needed to accomplish the placement 2,3 . 
However, full state feedback is not feasible for many control applications such as control of 
flexible structures. The dynamics of flexible structures is described by infinite dimensional 
equations and is usually approximated by large second-order discrete models; therefore, the 
measurement of the full state is not feasible. Output feedback pole placement is the practical 
choice for such systems. However, much of the literature for output feedback eigensystem 
assignment deals with relatively low order systems, where the number of states is less than 
the sum of inputs and outputs 4,5,6 , and therefore is not applicable for problems like control 
of flexible structures. Conditions on the closed-loop eigenvectors for the existence of output 
feedback gain matrices to accomplish the desired eigenvalue placement are discussed in Ref. 
7. References 8 and 9 present techniques for synthesis of output feedback gains, which assign 
a self-conjugate subset of the closed-loop eigenvalues. However, a major drawback of using 
output feedback for eigensystem assignment of large order systems is that while assigning a 
self-conjugate set of eigenvalues, the feedback may cause other closed-loop poles of the system 
to become unstable 10 . Therefore, it is desirable to develop techniques for output feedback gain 
matrices which assign dominant modes of the system while ensuring that the overall system is 
not destabilized. This problem is addressed in the current paper for strictly passive linear time- 
invariant systems in descriptor form, which are used in the characterization of the dynamics 
of flexible space structures. 

Robust control of flexible space structures is difficult due to its dynamic characteristics, 
namely, numerous closely-spaced, low frequency modes with little inherent damping 11 . 
Moreover, the parameters associated with flexible structure models such as modal frequencies, 
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damping and mode shapes often cannot be identified accurately. Feedback control is further 
complicated by instabilities induced by observation and control spillovers 11 . However, with 
collocated and compatible actuators and rate sensors, the system dynamics are passive, 
irrespective of unmodelled dynamics and parametric uncertainty 12,13 ’ 14 . Moreover, by the 
passivity theorem 15 , a strictly passive controller guarantees robust stability of the closed-loop 
system. For constant gain feedback, strict passivity requires that the symmetric part of the 
constant gain matrix be positive definite. Note that symmetry of the gain matrix is not necessary 
for strict passivity. Positive definiteness of the symmetric part of the output feedback gain 
matrix is necessary and sufficient for strict passivity. Similarly, when the plant is strictly 
passive, a passive controller ensures robust closed-loop stability. Positive semidefiniteness of 
the symmetric part of the gain matrix ensures passivity of a constant gain feedback controller. 
Thus, a strictly passive system is guaranteed to be robustly stable with constant gain output 
feedback if the symmetric part of the gain matrix is positive semidefinite. 

Thus, the robust eigensystem assignment problem for strictly passive systems is to construct 
output feedback gain matrices whose symmetric part is positive semidefinite. Preliminary 
approaches to these problems are described in Refs. 16 and 17. Reference 17 presents a 
sequential procedure through which a pair of self-conjugate eigenvalues is assigned to desired 
locations at every step with the aid of generalized Schur decompositions and orthogonal state 
transformations. This procedure is repeated until up to m self-conjugate closed-loop eigenvalues 
have been assigned, where m is the number of control inputs to the system. At each step 
the output feedback gains are designed such that the previously placed eigenvalues remain 
unchanged. This approach reduces the problem of assigning one self-conjugate pair of closed- 
loop eigenvalues to constrained solution of a set of quadratic equations. In this paper, necessary 
and sufficient conditions for the existence of a solution to these quadratic equations are presented 
in this paper for two cases, namely, the nonsymmetric gain matrices and the symmetric gain 
matrices. Furthermore, a minimax approach is described for obtaining parameters which satisfy 
the conditions for pole placement. Once parameters satisfying these conditions for pole placement 
have been determined, the approach to evaluate output feedback gain matrices follows. These 
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results are first developed for assigning dominant poles of passive LTI systems in descriptor 
form, in particular, for damping enhancement of flexible space structures with rate feedback. 
Then, these techniques are extended to pole placement with both position and rate feedback for 
second-order models of mechanical systems, while ensuring stability of the overall closed-loop 
system. The conditions guaranteeing stability in this case are that the symmetric part of the 
rate gain matrix is positive semidefinite, and the position gain matrix is symmetric and positive 
semide finite. Corresponding quadratic equations and conditions for existence of a solution for 
eigenvalue assignment using both rate and position feedback are presented. Finally, a numerical 
example of damping enhancement in the dominant modes of a structural testbed at NASA 
Langley (CEM phase II) is presented to demonstrate the application of this approach. 

Passive LTI Systems 

Systems that do not generate energy internally are referred to as passive systems. 
Mathematically, input-output systems are defined to be passive if (y,u) > I, where (•,•) 
denotes the inner product of the input-output vector spaces, u is the input to the system, y is 
the output of the system, and I is a parameter corresponding to the initial conditions of the 
system 15 . First, consider a constant gain system, y = Gu , where G is the constant gain matrix. 
Note that y T u = u T Gu = i/ T (sym{(S'}) i/, where sym{ G } denotes the symmetric part of G, 
defined as sym { G } = ^ ( G + G T ). Now, (//. u) > 0 (for passivity with / 3 = 0) if and only if 
sym {(',) > 0, where “> 0” denotes positive semidefiniteness of the matrix on the left-hand 
side of the inequality. Thus, it follows directly from the definition that positive semidefiniteness 
of the symmetric part of the gain matrix, G, is equivalent to passivity of a constant gain system. 

Next, consider passivity of LTI systems in descriptor form. Such LTI systems are expressed 
as 

Ez = Az + Bn, y = Gz + Du (1) 

where ~ is the system state, y is the system output, u is the system input, and the matrices 
(E, A , B, ( ' . D ) describe the dynamics of this system. For regular descriptor systems, the transfer 
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function matrix of the LTI system is G(s) = C(sE — A) -1 B + D. Passivity of an LTI system can 
equivalently be characterized as G(s) + G*(*) > 0 for all Re{s} > 0 15 . Extending the Kalman- 
Yakubovitch lemma for passive LTI systems (or equivalently, positive real transfer functions) to 
descriptor form, it follows that if there exist matrices, P = P T > 0, L and W, which satisfy 

a t pe + e t pa= -l t l 

E t PB-G t = -L t W (2) 

I) + I) 1 = W T W 

then the transfer function G(s) is positive real. Note that for E = /, the conditions in Eq. (2) 
reduce to the necessary and sufficient conditions of the positive realness lemma 12,13 . 

Dynamics of a flexible space structure with rate output provides an example of passive 
LTI systems. The dynamics of flexible space structures are typically expressed in second order 
form as 


Mx + Dx + Kx = Fu 


( 3 ) 


where x denotes an n x 1 vector of displacements, it denotes an m x 1 vector of control inputs, 
M, I). I\ are mass, damping and stiffness matrices of the system, respectively, and F is an 
n x m input influence matrix describing the actuator force distribution. The mass matrix, M, is 
a symmetric, positive definite matrix for physical systems, and the damping and stiffness matrices 
are symmetric, positive semide finite matrices. Using compatible rate sensors, collocated with 
the actuators, the m x 1 vector of rate output of the system is given by y = i' ] x. Compatibility 
of sensors and actuators implies that angular rate sensors are used with torquers, and linear 
velocity sensors are used with force or thrust actuators. These conditions ensure that the output 
influence matrix is the transpose of the input influence matrix. The dynamic equations for 
second order models can be rewritten as 
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that is, the descriptor form 



( 4 ) 


Ez = Az + Bu , y = Gz 


( 5 ) 
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where : = 'j. , E = q ^ , A = , B = ^ and C = B T . Note that, 

the conditions for positive realness in Eq. (2) above are satisfied with P = 

L = \ 0 x/2 A i) ] and W = 0, where I./) is a Cholesky factorization of D, that is, D = L L d . 

Note that the mass matrix is nonsingular, and the stiffness matrix must be nonsingular for 
complete observability with rate feedback. Thus, the matrices above satisfy the conditions of 
the lemma for passivity of descriptor systems. Note that symmetry and positive definiteness 
properties of the matrices M, I\ and D follow from the physics of flexible structure dynamics and 
do not depend on specific values of the mass, stiffness and damping parameters. Thus, dynamics 
of flexible structures with rate output are passive irrespective of the parametric uncertainties 
associated with the mass, stiffness and damping matrices of the system. Furthermore, truncation 
of the system dynamics to a smaller set of modes for reduced-order design models does not 
affect these passivity properties. 

The control law for constant gain output feedback is u = —Gy, where G is a constant m x m 
output feedback gain matrix. The closed-loop system dynamics with output feedback are given 
by E z = [A — BGB t )z. If the system is strictly passive, stability of the closed-loop system is 
guaranteed by the passivity theorem 15 for any output feedback gain matrix whose symmetric part 
is positive semidefinite. The main problem addressed in this paper is how to select a constant 
output gain matrix, whose symmetric part is positive semidefinite, for assigning dominant poles 
of a system to enhance closed-loop system response while ensuring over all closed-loop system 
stability. The sequential procedure for assigning desired closed-loop eigenpairs described in this 
paper assumes that the open-loop LTI system is strictly passive. Flexible space structures with 
zero frequency rigid body modes, or modes with no damping, do not satisfy this condition. To 
deal with such systems, static dissipative position and rate feedback controllers may be designed 
first to ensure strict passivity of an intermediate closed-loop system 12 , and the sequential 
approach may be applied to this system for eigensystem assignment. Furthermore, subsequent 
discussions assume that a low order control design model for plant dynamics, as opposed to 
the full order finite element matrices, is used in the control design process. The reduction may 
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be achieved through a number of model reduction techniques, such as modal truncation and 
modal cost analysis. 

Sequential Procedure 

The task of assigning dominant poles of strictly passive LTI systems with output feedback 
gain matrices whose symmetric part is positive semidefinite is accomplished in a sequential 
manner. In each step of the sequential procedure, one self-conjugate pair of closed-loop 
eigenvalues is assigned to desired values while making sure that the previously assigned closed- 
loop eigenvalues are not disturbed. The approach is an extension of the ideas in Refs. 9 and 
20 to descriptor systems with output feedback. Generalized ordered real Schur transformations 
of the system matrices, E and A, are employed in each step of the sequential procedure, using 
the following result from Ref. 21. Given two real matrices E and A, there exist orthogonal 
transformations L and U such that the matrices E = L T EU and A = L T AU are in ordered real 
Schur form. This transformation is known as the generalized real Schur transformation. These 
transformations are used to move previously assigned eigenvalue pairs to the top left block of the 
pair (^E, A^j , and the structure of the new gain matrix is prescribed such that it only affects the 
eigenvalues in the lower bottom portion of the system matrices. At each step, after computing 
the gain matrix that assigns a pair of desired closed-loop eigenvalues, the intermediate closed- 
loop matrix is transformed to a generalized Schur form with all previously assigned eigenpairs 
in the top left block of the updated system matrix. The overall output feedback gain matrix 
is constructed by accumulating the gains from each step. Since the symmetric part of each of 
the component matrices is positive semidefinite, the overall gain matrix also satisfies this robust 
stability condition. This process can be continued until m closed-loop eigenvalues have been 
assigned to the desired locations. 

The sequential procedure is conceptually straightforward and easy to implement in software 
since updated variables in the sequential process overwrite previous variables. However, 
description of the procedure below appears to be somewhat complicated due to notational 
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details in expressing the algorithm. For initialization of the procedure, denote the original state 
coordinates as | = ~, the system matrices in the original state coordinates as E\ = E, A\ = A, 
B\ = B, C i = ( \ and, initially, the output gain matrix is G 0 = 0. At the k th step of the 
sequential procedure, (k — 1) self-conjugate pairs of the closed-loop eigenvalues have been 
placed in the previous steps. A gain matrix is computed in this step to place the k th pair 
of the desired closed-loop eigenvalues while making sure that the (A — 1) self-conjugate pairs 
assigned in the previous steps remain unchanged. 

First, orthogonal matrices A /, and ( j are selected such that A’/ = l. [ E^l j and A /, = 
kj. Ai,l /, are in ordered real Schur form, and the {k — 1) self-conjugate eigenvalue pairs assigned 
in the previous steps are in the top left block of Applying these transformations to 

the system equations leads to /A/, = l. j. /A/ and ( '/, = ( '/, ( ), . If A/, is the gain matrix at the 
k th step, the closed-loop equations would be of the form /A/, A/, = ^/l/ — BjCi/ '/.j-/., where 
= / j/ :/, is the transformed state. The goal is to prescribe the structure for the gain matrix 
A'*; such that previously placed eigenvalues are not disturbed by this gain matrix. 

The matrices E^.A^ in ordered real Schur form are partitioned as follows: 


Ek = 


Em 

Ekl2 

, Aft = 

aU-ii 

Ak 12 

0 

Ek22 - 


0 

A-k > > . 


( 6 ) 


where E^n, A^ n are 2( k — 1) x 2( k — 1) matrices, E^, Akr 2 are 2( k — 1) x 2 ( ?? — k + 1) 
matrices, and A’/ _^. .1 /, _)_) are 2(?? — A + 1) x 2(?? — A + 1) matrices. Note that the orthogonal 
matrices A /, and / j have been selected such that the previously placed closed-loop eigenvalues 
are the eigenvalues of /A/. 1 1 . and the remaining eigenvalues of the system are 

eigenvalues of (e ^ 22 , .1 ) . Corresponding partions of /A/, and A j are 


Bk 


Bjk\ 

.Bk 2 . 


= [A,! C, 2 ] 


( 7 ) 


where AA fcl is a 2(A — 1) x ??? matrix, B^ is a 2(?? — A + 1) x m matrix, C kl is a m x 2(A — 1) 
matrix, and ( '/, 2 is a ??? x 2(?? — A + 1) matrix. 

Let Aj, denote a matrix whose columns form an orthogonal basis for the left null space of 
( '/, | , that is, Aj is a matrix with orthogonal columns such that A j/Y j | = 0. Let ??? denote the 



left nullity of C k \, so that A /, is an m x ???.*. matrix. If the gain matrix Gk is constructed as 
Gk = -V//// Aj/, where A/, may be an arbitrary ???.*. x ??? ^ matrix, then output feedback with A'/, 
as the gain matrix will not affect the (k — 1) eigenvalue pairs assigned in the previous steps. To 
see this, consider the expanded form of the following closed-loop matrix, 


Ak - BkG k C k 


Am Ak\2 — BkiNk^k^kCh-2 
. 0 Ak22 ~ Bk2NkGkN^Ck2 . 


(8) 


The first column of submatrices is not affected by the feedback since N£C kl = 0 by 


design. Therefore, eigenvalues of yEk, Ak — BkGkCk j are the union of the eigenvalues of 
(A/,i i • A k n) and the eigenvalues of (c k >>- A k22 ~ kC k^'J C k>^j , that is, the previously 


placed eigenvalues of ( A/ 1 1 . Akn ) remain unaffected by the output feedback. 


The matrix G k is selected to place the k th self-conjugate pair of closed-loop eigenvalues, 
(A*, A*). Denoting E k = E k 22 , A k = A k22 , B k = B k2 N k , and C k = A / ( the problem is 
reduced to selecting an ???.*. x nik output feedback gain matrix, G k for the descriptor system 
represented by the system matrices E k , Ak, Bk, and ( '/, , such that (A/, . A/, ) is an eigenvalue 
pair of the closed-loop system i^Ek, Ak — Ci/'i/ k^j ■ It must be assumed at this point that the 
controllable and observable subspaces for the system ^ A k , Bk, ( '/, j are not null. Note that, by 
construction, if the symmetric part of G k is positive semidefinite, then the symmetric part of 
C k will be positive semidefinite, and if G k is symmetric and positive semidefinite, then G k is 
symmetric and positive semidefinite. Thus, if A/, satisfies the criteria for robustness, so does 
A/,. The procedure for selecting G k is described in the next section. 


Finally, the gain matrix is updated as G k = G k -\ + Gk, and the system matrices are updated 
as E k+ 1 = E k , A k+ 1 = (A/, - Ih/nJ'k), B k+1 = B k , and C k+ i = C k , for the succeeding 
steps of the sequential procedure. Note that ???.*., the dimension of the left null space of Cki 
decreases by 2 in each step of the sequential procedure, that is, m k = m-k - 1 — 2. 


The sequential procedure above is repeated until all the desired dominant modes have been 
shifted, up to a maximum of m/2 modes or complex-conjugate pairs, where m is the number 
of inputs. After m/2 iterations of the procedure described above, m/,, the dimension of the left 
null space of ( /, | will be zero, and there would be no freedom to assign additional eigenvalues. 
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Eigenpair Assignment 


This section describes the approach to select output feedback gains to assign one pair of 
complex conjugate eigenvalues, while ensuring that (1) the symmetric part of the gain matrix 
is positive semidefinite, or (2) the gain matrix is symmetric and positive semidefinite. For 
notational simplicity, the system matrices will be denoted as E, A , IL C (corresponding to 
. ,1 /, . /!/ . ( '/, of the previous section), the output feedback gain matrix will be denoted as 
G (instead of G /, of previous section), the number of inputs is denoted by m (in place of ////,) 
and the desired eigenvalue pair will be denoted (A, A) rather than (A k , A /, ) . Therefore, in the 
simpler notation, the problem is to select a matrix G, with a specified structure for robustness, 
such that (A, A) is a generalized eigenpair of the closed-loop system matrix, (E, A — BGG ). 

Let o be the closed-loop eigenvector corresponding to the eigenvalue A. The generalized 
eigenvalue problem becomes (A E — A + BGC)4> = 0. This closed-loop expression can be 
rewritten as 


[A E - A 


B 


= r 


= n 


(9) 


GCcb J L GCcb J 

It is obvious from Eq. (9), that the vector on the right hand side of the expression above must 
lie in the right null space of T. Let A be a matrix whose columns form an orthogonal basis 
for the null space of T, that is, TiV = 0. Note that though E,A and B are real matrices, T 
and N are complex matrices since the eigenvalue A is complex scalar. However, to ensure 
that the gain matrix is real the closed-loop eigenvector corresponding to the complex-conjugate 
eigenvalue is chosen to be the complex-conjugate of <£>, that is, o is chosen to be the eigenvector 
corresponding to A. The overbar in the expressions in this section refer to complex-conjugation 
of the elements of the corresponding vector (or matrix) only, as opposed to the Hermitian 
operator, which involves transposition and complex-conjugation. 

Since columns of N span the null space of T, it follows that 


</> 

= Na = 

- N\ - 

_ GC<t> _ 


No 


o 


( 10 ) 
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where a is an arbitrary vector of complex elements, and the matrices .Y| . A4 are formed by 
partitioning N compatibly with o and GC<f>. Therefore, <f> = N\a and GC<f) = i\4c% which 
leads to 


GCN\ot = N -2 a (11) 

The eigenassignment problem is now reduced to selecting a such that there exists a gain 
matrix, G, satisfying Eq. (11) whose symmetric part is positive semidefinite. With o being the 
eigenvector corresponding to A, real solutions for the gain matrix G can be obtained, and the 
equations can be written out to involve only real arithmetic operations as follows. 

For the eigenvalue, A, with closed-loop eigenvector, <f>, the matrix T = [\E — A \ B], 

and N is a matrix whose orthogonal columns span the null space of T. If the arbitrary coefficient 
vector is chosen to be a, the complex-conjugate of «, then it follows that 

GCNia = N 2 a ( 12 ) 


Eq. (11) and Eq. (12) can equivalently be rewritten as 

Re(o) 
Im(o) 


r;C[Re(iVi) -Im(iVi) 


Re(o) 

Im(o) 


= [Re(iV 2 ) -Im(iV 2 ) 


(13) 


and 


GC [ Im( iVi ) Re( iVi ) 


R.e(o) 

Im(o) 


= [Im(iV 2 ) Re(iV 2 ) 


Re(o) 

Im(o) 


(14) 


where Re(«) denotes real part of the argument, and 1 1 n (•) denotes imaginary part of the 
argument. In compact form these equations are written out as 

GlV lP = V lP 

gw 2P = V 2P 


(15) 


where P = [Re(o); Im(o)], ITA = C[Re(iVi) -Im(iVi)], Vj = [Re(iV 2 ) -Im(iV 2 )], 
W 2 = f '[[in(Aj) Hc'( Aj ) ], and 14 = [lin( A 2 ) H<'( A 2 ) ]. Note that Eq. (15) is a system 
of quadratic equations in the unknown variables, namely, the elements of the gain matrix, G, 
and the coefficient vector, />. Further, the elements of G are constrained such that sym(G') > 0 . 
Solution of Eq. (15) provides the desired pole placement with dissipative gain matrices. The 
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following proposition gives the necessary and sufficient conditions to be satisfied by a coefficient 
vector p for the existence of a solution matrix G of desired structure 22 . 


Proposition 1. There exists a matrix G, whose symmetric part is positive semidefinite, that 
satisfies Eq. (15), if and only if, there exists a vector p which satisfies 


P T \fw lP > n 
p T \fw 2 p > n 


(16) 


P T \f w 1P ^ (p T v?w 2 p) - \(fv? w 2 p + p T \fw lp y > n 


Furthermore, there exists a symmetric, positive semidefinite matrix G which solves Eq. (15) if 
and only if there exists a vector p which satisfies Eq. (16) and p T \fw 2 p = p T vTw 1 p. 


Proof: Denote y\ = l) />. y> = V 2 p, x \ = W\p and x 2 = W 2 p. Then, Eq. (15) leads to 
G\x | x 2 ] = [; f/i y 2 ] , and the conditions in Eq. (16) become yfx i > 0, y 2 x 2 > 0 and 
{yTzi){y 2 Z 2 ) ~ \{yTx2 + y 2 x iy > 6- Therefore, 


[.ri 


x 2 ] T G[x i x 2 


[.ri x 2 ] T [y 1 y 2 ] 

r T T 1 

y i y-2 

T T 
|_.r 2 y i x 2 y 2 \ 


(17) 


Adding this relation and its transpose yields 


T 

[•?T ^2] [sym( (S')][.ri .r 2 ] = 


T 

y 1 ^t 


0.5 (yf x 2 + y 2 xi 


(18) 


Now if sym(G') > 0, then the equation above implies, 

VlX x 0.5 ( yj x 2 + yT x 1 ) 

0.5 (yf x 2 + yT ;ri ) y 2 x 2 


> 0 


(19) 


Noting that determinants of principal minors of a positive semidefinite matrix must be 
nonnegative leads to conditions in Eq. (16). Thus, if there exists a matrix G, which satisfies Eq. 
(15), and its symmetric part is positive semidefinite, then the vector p must satisfy conditions 
in Eq. (16). 


Conversely, it has to be shown that if a vector p satisfies the conditions in Eq. (16), then 
there exists a matrix G which satisfies Eq. (15), whose symmetric part is positive semidefinite. 
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This is proved by constructing such a matrix. Let X = [,ri x> ] and Y = [y\ y> ] , then Eq. 
(15) is rewritten as GX = TA Let Q be an orthogonal matrix, such that 

Q t Y = ] ( l (20) 

where l) is a nonsingular 2x2 matrix. The matrix Q can be obtained by QR factorization 
of 1'. Now, define A’ | . AA as follows 

| 1 ]=Q t X (21) 

_ A 2 . 

where A’| is a 2 x 2 matrix, and AA is a (m — 2) x 2 matrix. Now A’| is nonsingular if x \ 
and x -2 are linearly independent (otherwise, Eq. 15 is solved trivially). Defining G\\ = V) A’ ^ 1 , 
it can be seen that 



0 


0 




( 22 ) 


Therefore, it follows that the matrix G defined as G = Q 1 jj () l satisfies GX = Y. 

Next it is shown that sym(G') > 0. From Eq. (19) it follows that conditions in Eq. (16) 
imply sym (Y T X) > 0. Next, since Y r X = ((.)' ) ')' ((.)' A ) = ) \' A , , sym(F T X) > 0 leads 
to sym ^V'i^A’i j > 0. Noting that sym(Yi T A'i) = Xf Jsym^&' n ^ J > 0, sym^G'n^ > 0 

follows from the nonsingularity of A | . Finally, by construction, sym(G') > 0 if sym^&'ii^ > 0. 

Furthermore, if G is symmetric, that is, G = sym(G'), then yfx 2 = x^G T x> = x^G T x 1 = 
y^x 1 , so that the additional condition in the proposition is satisfied. On the other hand, 
the additional constraint ensures that Y T X = V'/ A’i is symmetric. Since Y\ = G \ \ X \ , 
Y?X 1 = Xf l'i , and A', is nonsingluar, it follows that G \ \ is symmetric. Finally, again 
by construction, G is symmetric when G \ \ is symmetric. ■ 

With this proposition, the problem of eigenpair assignment with dissipative output feedback 
is reduced to selection of a coefficient vector p which satisfies the conditions in Eq. (16). A 
minimax approach to determining a coefficient vector p which satisfies the conditions in Eq. 
(16) is described next. 
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Note that if there exists a coefficient vector p which satisfies 


P T \fw lP - 

p T \fW 2 p- 

then it also satisfies the conditions in Eq. (16). The first two inequalities of Eq. (16) are 
obviously satisfied when the inequalities in Eq. (23) are satisfied, and multiplying the two 
inequalities in Eq. (23) shows that the third condition of Eq. (16) is also satisfied. Further, the 
inequalities in Eq. (23) can equivalently be written as four inequalities, which are quadratics 
in the coefficient vector, p, as follows 

fi(p) = /{if fEi + 1 -(V?W 2 + Vf Wi) jp > 0 

hip) = P T Ufw 1 - \(v?w 2 + vhwMp > 0 

; (24) 

hip) = P T \vhw 2 + - (i'i T it 3 + \fw x ) > p > n 

hip) = p T {vf w 2 - l -(y?w 2 + vf Wi) |p > 0 

A feasible coefficient vector, p, can now be obtained maximizing the minimum of the four 
quadratics, fi{p), | = 1, ...,4. This problem becomes a standard minimax problem by reversing 
the sign of the quadratics, that is, minimizing the maximum of —fi{p), i = This 

problem is transformed to a constrained minimization problem by introducing a scalar variable, 
7], as follows 22 , 


\{p T v^w 2 p + p T \fw lP ) 

^ (p T Vi W 2 p + p T Vh Wipfj 


> n 

> n 


( 23 ) 


min 7] such that flip) + 7] > 0, i = 1, ...,4 (25) 

p,v 

Standard nonlinear programming techniques may be used for this constrained minimization. 
Analytic gradients of f , ( /> ) are readily available, since the gradient of any quadratic, f{p) = 
p T Qp, is given by = (Q + Q T )p. The optimization of Eq. (25) is very well-behaved 

because the functions, fi{p), i = 1, ...,4, are quadratic in p, and analytic gradients are linear in 
p. A nonpositive value of // in Eq. (25) provides a feasible coefficient vector, p, that satisfies 
the conditions in Eq. (16) for existence of a gain matrix whose symmetric part is positive 
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semidefinite. The convergence of the nonlinear optimization in Eq. (25) is not an issue, since 
the search can be terminated once a desirable target value of // (which must be negative) has 
been attained. Furthermore, if a symmetric, positive semidefinite gain matrix is desired, the 
additional equality constraint required for symmetry (in Proposition 1) can also be included in 
the optimization of Eq. (25). Since this equality again involves a quadratic in p, its analytic 
gradients are linear and readily available. 

Experience in application of the minimax approach presented above has shown that this 
technique is very effective in obtaining a feasible coefficient vector, p, which satisfies the 
conditions in Eq. (16). However, a number of other approaches have also been explored 
for this problem 16,17 . Another useful approach is to characterize a subspace for the vector p, 
which satisfies the first two conditions in Eq. (16). This can be accomplished by examining 
the eigenvectors corresponding the positive eigenvalues of the symmetric parts of the matrices 
l | / ff'i and l'/ff,. Then the function in the third inequality of Eq. (16) is maximized with 
p being restricted to this subspace, until a positive value for that function is reached. This 
approach has also proved successful for a number of problems, though the minimax approach 
presented earlier has been successful in obtaining a feasible coefficient vector, p, for a larger 
class of problems. 

Once such a coefficient vector has been determined, the procedure given in the proof 
of Proposition 1 can be followed to determine the gain matrix G, with positive semidefinite 
symmetric part, that assigns the desired closed-loop complex-conjugate eigenvalues. 


Position and Rate Feedback 

Pole placement for second order models, such as those for dynamics of flexible space 
structures, with both position and rate feedback is considered in this section. Position feedback 
gives additional freedom to place closed-loop eigenvalues, which may not be feasible with rate 
feedback alone. The technique for eigenpair placement with rate and position feedback parallels 
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that for rate feedback alone, the difference being that the development is more involved since 
expressions for both rate and position feedback gain matrices are considered simultaneously. 

Consider a second-order model, as shown in Eq. (3), again. Position sensors are now 
assumed to be collocated with the rate sensors and actuators to provide the position output. 
Thus, an m x 1 vector of position output, y p , is given by y p = C p z , where C p = [F t o] . As 
before, the rate output is //, = ( : , where C r = [o F T ]. 

The eigensystem assignment problem is to determine an m x m symmetric, positive 
semidefinite position gain matrix, G p , and an m x m rate gain matrix, G r , whose symmetric part 
is positive semidefinite, such that the control law u = —G p y p — G r y r assigns dominant poles 
of the system to desired locations. In other words, to determine gain matrices G p and G r such 
that certain eigenvalues of {E, A — BG p G p — BG r C r ) are at desired locations. Again, it must 
be ensured that other closed-loop eigenvalues do not become unstable during output feedback 
eigensystem assignment with position and rate feedback. Using Lyapunov function arguments, 
it can be shown that if G p is a symmetric, positive semidefinite matrix, and the symmetric part 
of G r is positive semidefinite, then the closed-loop system is stable. 

A sequential methodology is used for eigensystem assignment with position and rate 
feedback, extending the approach for rate feedback in previous sections. At each step of 
the process, one self-conjugate eigenpair is assigned, while employing a structure for the 
gain matrices such that previously placed eigenvalues are not disturbed. For initializing the 
sequential procedure, set E\ = E, A\ = A, B Pt i = B, B r> \ = B, C p> \ = C p , C r , 1 = C r , 
the original coordinates as ~i = ~, and, the output gain matrices as G 1P 0 = U,.o = 0. At 
the k th step of the sequential procedure, orthogonal matrices L k and I j are selected such 
that Ek = l- j. lyj j and .1 /, = /.J.l/V j are in ordered real Schur form, and the {k — 1) self- 
conjugate eigenvalue pairs assigned in the previous steps are in the top left block of E/ . Aj^j . 
Applying these transformations to the system matrices leads to B p y. = l. j. H p /^ B r y. = kj. /i, / , 
C p ,k = C P ,kUk and C r ,k = Cr.kUk ■ Partition C pk and C r k as follows: C pk = C pM C Pjk2 
and C r k = C r kl C rk2 , where G pkl , C r k i are m x 2(k — 1) matrices. System matrices 

A k . Il pk and B rk are also partitioned conformally. Let the columns of N pk form an 
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orthogonal basis for left null space of C pk i , and let the columns of .Y, / form an orthogonal basis 
for left null space of C Tj ki . Form the rate and position gain matrices for the k th step as follows: 
G pk = Np^Gp^Np/., where G pk is an arbitrary matrix to be defined; and G r , k — N Tj k G r , k , 
where G r k is an arbitrary matrix to be defined. Then, the matrices G p k and G r k do not affect 
the previously placed eigenvalues. To see this, consider the expanded form of the closed-loop 
system matrix, (/\ k — B pk G pk G pk — B rk G rk C rk \, shown below. 


YU-ii 

0 


Yifci2 B pk \N pk G pk N p ^ k G ]) : k2 B rk \N rk G rk N r , k G rk 2 
A k 22 B pk 2Np k Gp k N^j,C'p k 2 B rk 2N rk G rk N r j,C rk 2 


( 26 ) 


The first column of submatrices remains unaffected by the feedback because 
N pk C p , k 1 = 0 and ' rk | = 0 . Thus the eigenvalues of the closed-loop system 

with position and rate feedback of this form are the union of the eigenvalues of 

(^E k iii Ylm), which are the previously placed eigenvalues, and the eigenvalues of 

( K B k 22,A k 22—Bp k 2Np k Gp k Nj k Cp k 2—B rk 2N rk G rk N^ k C rk 2^. Denoting E k = E k 2 2, 
A k — Yl^22, B pk — Bp k 2Np k , B rk — B rk 2 N rk , C pk — Gk^ k Cp k 2, and C rk — Gf rk f rk -2, 
the problem is to select matrices G pk and G r k that satisfy the closed-loop stability 
constraints and that (^E k , A k — B pk G pk G pk — B rk G rk C rk ^J has one self-conjugate pair 
at the desired values. The gain matrices are updated as G pk = G pk _i + G pk and 

G , k = G ,.k- 1 + G l k , and the system matrices are updated for the succeeding steps as: 

B k ~\- 1 — E k , A k _ |_^ — (yA k Bp k Gp k Cp k B rk G rk C rk ^j, B pk _^i — B pk , B rk _ |_i — B r k , 
G p k+ i = O p k and ( ' l k . | = G r k . This procedure may be repeated up to m /2 times, where 
m is the number of inputs to the system. 

For notational simplicity, consider synthesis of position and rate gain matrices G p . G , for 
eigenplacement of the system given by the matrices E, A, B p , B r ,C p ,C r . In other words, the 
problem has been reduced to selecting a symmetric, positive semidefinite position gain matrix, 
G p , and a rate gain matrix, G r , whose symmetric part is positive semidefinite, such that (A, A) 
is an eigenpair of the generalized eigenvalue problem (E, A — B p G p Op — B r G r G r ). Let <f> be 
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the closed-loop eigenvector such that (A E — A + B P G P C P + B r G r C r )(f> = 0, or equivalently, 


[A E- A 

1 Bp 

1 B r ] 

G p Cp(f) 

= r 

GpC p (j) 




G r C r <f>_ 


G r C r <f>. 


(27) 


Let N be a matrix whose columns form the right null space of T. Then, solution to Eq. (27) 
may be written as 




' Ni ' 

GpC p <t> 

= 

N p 

G r C r <f>_ 


N r 


(28) 


where A) . N p and N r are conformal partitions of the matrix A , and a is an arbitrary coefficient 
vector. Eq. (28) leads to 


GpCpNia = N p a 
G r C r N ia = N r a 


(29) 


Similar treatment for the complex-conjugate eigenvalue, A, and a corresponding eigenvector, 
o, leads to 


G P C P N iq = N p a 

(30) 

GrCrNia = N r a 

Combining these equations to restrict all operations to real arithmetic results in 

G r W r ip = V r ip GpWpip = Vpip 

(31) 

G r W r -2P = V r 2P G p W p 2P = V p 2 p 

where p = [Re(o); Im(o)]; W rl = C, llci Aj ! -Im(iVi)]; V rl = [Re(iV r ) -Im(iV r )]; 
W r2 = CV[Im(iVi) Re(iVi)]; fy 2 = [Im(7V r ) Re(iV r )]; W pl = ( ' p \iv\ Aj ! -Im(iVi)]; 
Vpi = [Re(iVp) -lm(Np)]; W p2 = C p [Im(iVi) Re(iVi)]; and V p2 = [Im(iV p ) Re(iV p )]. 

Thus, the problem of eigensystem placement with rate and position feedback is reduced to 
selecting a coefficient vector p, a symmetric, positive semide finite matrix, G p , and a matrix, 
G r , whose symmetric part is positive semidefinite which satisfy Eq. (31). The conditions for 
existence of a solution to this system of quadratic equations is given in the following proposition. 
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Proposition 2. There exists a symmetric, positive semidefinite matrix, G p , and a matrix, 
whose symmetric part is positive semidefinite, which satisfy Eq. (31), if and only if, there 
exists a vector p that satisfies 

P T v£w pl p > n p T v^w rl p > o 

p t v 1 J 2 w P 2p > n p T v^ 2 w r2 p > o 

(p T V^W p1P ) (p T V^W p2P ) - + P T V^W plP y > n (32) 

{fv^w rl pj (p T v£w r 2 p) - l -(p T v^w r2 p + p T v£w rlp y > o 
P T v£w p 2P = p T v£w pl p 

The proof of this proposition follows that of Proposition 1 for rate feedback. 

Again, a minimax approach can be employed to determine a feasible coefficient vector, p, 
which would satisfy the conditions of Eq. (32). It can be easily verified that a vector p satisfies 
these conditions if it satisfies the following eight quadratic inequalities, 

Mp) = + \(v^w p2 + v£w pl )} P > n 

hip) = p t {v£w p1 - \(y^w p2 + v^w pl )^p > n 

Mp) = p t {v£w P 2 + \{v£w p2 + vj 2 w pl )}p > n 

Mp) = p t {v£w p2 - \(v£w p2 + v£w pl )\ P > n 

r i i (33) 

Mp) = p T W£w r i + -(y£w r2 + Vr 2 w rl ) L > o 

Mp) = w rl - ^(v£w r2 + v£w r i ) |p > o 

Mp) = + l -(v^w r2 + v^w r i)^v > o 

Mp) = p T {v£w r2 - l -{v^w r2 + v^ 2 w r i ) }p > o 

and the equality constraint ensuring symmetry of the position gain matrix. A vector p satisfying 
these conditions can be determined by maximizing the minimum of fi{p), i = 1, ..., 8 . This 
is accomplished by performing a constrained minimization, similar to that of Eq. (25), with 
the additional equality constraint for symmetry of the position gain matrix. This minimax 
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approach for determining feasible coefficient vectors p satisfying the conditions of Eq. (32) 
has been very successful in practical application, though a number of other techniques have 
been attempted 16,17 . Once a coefficient vector satisfying the constraints in Eq. (32) has been 
determined, the procedure to construct the desired matrices also follows the same approach as 
for rate feedback. 


Numerical Example 

The approach for eigensystem assignment with dissipative gains has been applied for 
synthesis of controllers to enhance the damping of the phase 2 CSI Evolutionary Model (CEM), 
a testbed for control of flexible space structures at NASA Langley. Damping enhancement 
control for this structure represents a typical problem in the control of flexible space structures, 
since this structure possesses numerous, closely spaced, low-frequency modes, with uncertain 
modal parameters. Damping enhancement is required since inherent damping of the structure is 
very low. However, general output feedback pole placement techniques for increasing damping 
in the dominant low frequency modes have a tendency to destabilize the remaining modes of 
the system. When collocated rate sensors and actuators are used, the dynamics of the structure 
are passive, and the technique described in this paper can be applied for damping enhancement 
of the dominant modes, while ensuring stability of the overall closed-loop system. 

The phase 2 CEM structure consists of a 62-bay central truss (each bay is 10 inches long), 
along with two horizontal booms for suspension, a vertical laser, and a vertical reflector tower, 
as shown in Fig. 1. This structure has 10 modes with frequencies up to about 5 Hz., and 
95 modes with frequencies under 60 Hz. The first six modes are rigid body modes, due to 
suspension of the structure from the laboratory ceiling, that have frequencies up to about 0.3 Hz. 
Eight control stations housing collocated and compatible sensors and actuators are located at the 
bays shown in Fig. 1. Air thrusters providing linear forces are available at these locations, along 
the directions shown in Fig. 1, providing an S x 1 control input vector, u. Linear velocities are 
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assumed to be available at these locations, along the same directions, giving a 8 x 1 measurement 
vector, //. The problem is to determine an 8 x 8 rate gain matrix, G, whose symmetric part 
is positive semidefinite, such that the control law u = —Gy places a subset of the closed-loop 
eigenvalues at desired locations without destabilizing the remaining modes of the structure. 

A control design model of the first 10 modes of the structure is used for this numerical 
example. Low inherent damping ratio of 0. 1 percent has been assumed for the open-loop system. 
The open-loop eigenvalues along with damping and frequencies are shown in Table 1. Modal 
displacements at the eight sensor/ actuator locations are shown in Table 2. Since m = 8 control 
inputs are available, it is possible to assign desired amount of damping to up to 4 modes. 

The first objective was to increase damping ratio in the first mode to 10 % that is, the first 
desired pair of closed-loop eigenvalues was Ai^ = —0.0818 ± 0.8139 /. Proceeding as described 
above, the system of quadratic equations in Eq. (15) was constructed for this eigenassignment. 
Then, a coefficient vector, p, which satisfies the conditions in Eq. (16) was determined as 

p= [1.0, 1.0, 1.0, -1.0, -1.0, 1.0, -1.0, 1.0, 1.0, 

-0.2757, 0.8621, -1.0, 1.0, -1.0, 1.0, -1.0] 


using the minimax approach. Upper bounds on components of the coefficient vector, p, were 
set to 1.0, and lower bounds were set to —1.0, for the constrained optimization problem in 
determining this coefficient vector. Proceeding as described in the proof of Proposition 1, a 
gain matrix G \ for this eigenplacement was computed. Nonzero eigenvalues of the symmetric 
part of G\ are 0.4129 and 0.8035, which demonstrates that the symmetric part of G \ is positive 
semidefinite. Closed-loop eigenvalues with this gain matrix are shown in Table 3. It should be 
observed that the gain matrix, Gi, successfully increased damping ratio of the first mode to 10 
% without destabilizing other modes. In fact, damping is increased in most other modes due 
to the dissipative nature of the controller. 

For the second step, damping ratio in the second mode of the system was to be 
increased to 10 % resulting in the second pair of desired complex-conjugate eigenvalues to 
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be A 3 4 = —0.0830 + 0.8259j. A coefficient vector that satisfies the necessary and sufficient 
conditions for eigenplacement in Eq. (16) was determined as 

p= [1.0, 1.0, 1.0, 0.7088, 1.0, -1.0, -1.0, -1.0, -1.0, 

-0.5250, -1.0, 1.0] 


Note that the dimension of the coefficient vector, p, for the second step has decreased by 4 as 
compared to the first step, since the dimension of the left null space of C 12 for the second step 
has decreased by 2, that is, m 2 = 6 . The cumulative gain matrix Go for placing two eigenpairs 
has nonzero eigenvalues of the symmetric part at 0.2090, 0.6899, 0.9358 and 1.8008, which 
demonstrates the stability constraint on the structure of gain matrices is satisfied. The closed- 
loop eigenvalues with the gain matrix, Go, are presented in Table 4, showing that both pair of 
complex conjugate eigenvalues are successfully placed without destabilizing any other modes. 

In the third step, damping ratio of the third mode was increased to 10 % which leads to 
the third pair of complex-conjugate eigenvalues as As^ = —0.0857 + 0.8522). Repeating the 
operations for the third step, a coefficient vector which satisfies the conditions in Eq. (16) was 

p= [-1.0, 1.0, -1.0, -0.9664, -1.0, 1.0, -1.0, 0.5762] 

which led to the cumulative gain matrix of G.\ for the third step with nonzero eigenvalues 
of the symmetric part at 0.0961 , 0.2843, 0.7060, 0.9694, 1.9922 and 4.2001. Closed-loop 
eigenvalues with gain G.> are shown in Table 5. Again, damping ratio in the first three modes 
has been increased exactly to 10 % as desired without destabilizing any other modes. In fact, 
damping in most other modes has increased, with mode 6 being overdamped. Furthermore, 
modes 4, 5, 7 and 8 are also very well damped. 

With 8 control actuators, it is possible to increase damping in one other mode. Therefore, 
mode 9 was chosen for an increase in its damping ratio to 10 % that is, the final pair of desired 
closed-loop eigenvalues was chosen as A 7 ,x = — 1.8692 ± 18. 5982j. A feasible coefficient vector 
for this case was determined to be p = [—1.0, —1.0, —0.7275, 0.0497]. The final gain matrix, 
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G, that places four desired closed-loop eigenpairs, without destabilizing other higher frequency 


modes, was computed to be 



- 0.7471 

0.2028 

- 0.3887 

- 0.8128 

- 0.2086 

0 . 116.3 

- 0.3025 

- 0 . 9071 - 


0.3282 

1.1480 

0.2282 

- 1.4228 

0.9391 

0.1604 

- 1.1181 

- 0.1716 


- 0.5818 

- 0.4845 

1.1284 

1.5011 

- 0.1217 

- 0.4460 

0.0589 

0.2609 

G 4 — 

- 0 . 4.390 

- 0 . 879.3 

- 0.3065 

3.1249 

- 1.3062 

0.4245 

1 . 19.33 

0.2856 

0.4710 

1.2129 

- 0.7990 

- 1.4584 

4.1617 

0.4467 

-. 3.9649 

0.1294 


0.4000 

0.1577 

- 0 . 5.367 

- 0.8686 

- 0 . 106.3 

0 . 381.3 

0.1086 

- 0 . 315.3 


- 0.7939 

- 1.3759 

1.1685 

1.7430 

- 4 . 158.3 

- 0.6551 

4.5028 

0.2142 


_— 0.8909 

0.0902 

0.3491 

0.3117 

0.3107 

- 0.0493 

0.4911 

1.2872 _ 


Eigenvalues of the symmetric part of the overall gain matrix, G, are 

0.0043, 0.0969, 0.1242, 0.5703, 1.0164, 1.9239, 3.0166 and 9.72S9, that is, sym(G) > 0. 
Closed-loop eigenvalues with the overall gain matrix is shown in Table 6, showing that all four 
closed-loop eigenvalue pairs have been placed at their desired locations, without destabilizing 
any of the higher frequency modes. In fact, damping has tended to increase in most other 
modes. This example clearly demonstrates the effectiveness of the proposed dissipative 
controller design technique for robust eigensystem assignment. 

Results of damping enhancement control with both position and rate feedback are presented 
next. First, damping ratio of mode 1 was increased to 10 %. Following a procedure outlined in 
the previous section, a symmetric, positive semidefinite position gain matrix, were determined to 
assign the closed-loop pair. Closed-loop eigenvalues with these rate and position gain matrices 
are shown in Table 7. Proceeding with the sequential procedure, the damping ratios of modes 
2, 3 and 7 were increased to 10 %. The eigensystem assignment procedure resulted in the 
following cumulative symmetric, positive definite position gain matrix, 


- 3.4841 

- 0.1125 

- 4 . 665.3 

0.2965 

- 0.2298 

. 3.1246 

1.1816 

- 1 . 9022 - 

- 0.1125 

1.9084 

- 0 . 58.32 

- 1.6722 

7.2152 

0.1647 

- 6.5656 

- 0.1680 

- 4 . 665.3 

- 0 . 58.32 

7.9072 

0.3844 

1 . 58.34 

- 4 . 2.352 

- 3.1930 

2.7889 

0.2965 

- 1.6722 

0.3844 

1.7415 

- 6.4787 

0.1705 

6.4519 

- 0.1206 

- 0.2298 

7.2152 

1 . 58.34 

- 6.4787 

46.6525 

- 0.2519 

- 45 . 6.357 

0.2129 

3.1246 

0.1647 

- 4 . 2.352 

0.1705 

- 0.2519 

3 . 34.34 

1.2556 

- 1.8501 

1.1816 

- 6.5656 

-. 3 . 19.30 

6.4519 

- 45 . 6.357 

1.2556 

46.9748 

- 1.4609 

_— 1.9022 

- 0.1680 

2.7889 

- 0.1206 

0.2129 

- 1.8501 

- 1.4609 

1.4274 _ 
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and a rate gain matrix, whose symmetric part is positive definite. 


- 0.6007 

- 0.7380 

- 0.6090 

0.7122 

0.4598 

0.4789 

0.4392 

0 . 0.355 

- 0.4611 

- 0.3645 

1.0371 

0.5755 

- 0.7028 

- 0.7606 

- 0.3794 

0 . 8 . 39.3 

0.9575 

- 0.3045 

- 0 . 42.36 

- 0 . 554.3 

0.2021 

- 0 . 7.357 

- 0.9152 

1.4169 

0 .. 359.3 

- 1.0655 

- 0.1207 

0.9327 

_— 1.2277 

1.3110 

0.0893 

- 1.6916 


- 0.4572 

0 . 6.321 

- 1.0302 

0.6469 

1.6212 

0.5455 

0.0844 

- 1.3040 

- 0.3726 

- 0 . 1.362 

0.9946 

0.3222 

- 0.6959 

- 1 . 4 . 3.34 

0.3812 

1 . 65.30 

. 3.5101 

0.5159 

- 3.4491 

- 0.6720 

- 0 . 096.3 

0.6254 

- 0.6292 

0.3972 

- 2 . 6.389 

0.1686 

3.3490 

- 0.4591 

0 . 22 . 3.3 

— 0.6664 

0.5661 

0.3660 


Table 8 shows the closed-loop system eigenvalues using the rate and position gain matrices 
obtained. It is observed that the desired closed-loop eigenvalues have been placed while 
maintaining stability of the other modes. Comparing the frequencies of the closed-loop poles, it 
is noted that some of the system frequencies have been shifted significantly from their open-loop 
values. Position gain brings about this shift in frequencies of the open-loop system, as may be 
anticipated. On the contrary, frequencies of the closed-loop system with rate feedback alone, in 
Tables 3-6, remain close to the open-loop values. This demonstrates one of the benefits of using 
both position and rate feedback, namely, significant shift in open-loop frequencies is possible 
if desired. Moreover, inclusion of the position feedback provides additional freedom beyond 
eigenvalue assignment such that improved controller designs are possible. 
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Concluding Remarks 


This paper presents a novel approach for robust eigensystem assignment using constant gain 
output feedback controllers with dissipativity constraints, that is, output feedback eigensystem 
with gain matrices whose symmetric parts are positive semidefinite. Robust eigensystem 
assignment in this context implies that stability of the closed-loop system is guaranteed, whereas 
in general eigensystem assignment using output feedback may destabilize the closed-loop system. 
This approach can be used for robust output eigensystem assignment for any passive LTI system 
in descriptor form. Specifically, it has been applied to second-order models of flexible structures 
with rate feedback. A sequential procedure is used to place one pair of complex conjugate closed- 
loop eigenvectors in each step of the procedure, while ensuring that the previously assigned 
eigenvalues are not disturbed using generalized Schur transformations. The problem of assigning 
one pair of complex-conjugate closed-loop eigenvalues with gain matrices of a specified structure 
is reduced to a constrained solution of a system of quadratic equations. Necessary and sufficient 
conditions for the existence of a solution to these equations are presented in terms of inequality 
constraints on a coefficient vector. A minimax approach is presented to determine coefficient 
vectors which satisfy these constraints. It is shown that the sequential procedure can be used 
to assign up to m self-conjugate closed-loop eigenvalues, where m is the number of control 
inputs. For second-order models of flexible space structures, the approach has been extended 
for robust eigensystem assignment with both position and rate feedback. A numerical example 
of damping enhancement for a large flexible structure has been presented to demonstrate the 
approach. This example clearly demonstrates that the proposed approach provides a practical 
technique to design controllers that enhance the performance of the open-loop system while 
ensuring that the overall system remains stable. 
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Open-loop Eigenvalues 

Damping Ratio 

Frequency 

(rad/sec) 

-0.0008 ± 0.8180 j 

0.0010 

0.8180 

-0.0008 + 0.8301 j 

0.0010 

0.8301 

-0.0009 + 0.8565 j 

0.0010 

0.8565 

-0.0011 + 1.1308 j 

0.0010 

1.1308 

-0.0011 + 1.1401 j 

0.0010 

1.1401 

-0.0019 + 1.9100 j 

0.0010 

1.9100 

-0.0107 + 10.7278 j 

0.0010 

10.7278 

-0.0149 + 14.9425 j 

0.0010 

14.9425 

-0.0187 + 18.6919 j 

0.0010 

18.6919 

-0.0341 + 34.0618 j 

0.0010 

34.0618 


Table 1. Open-loop eigenvalues 


Mode 

No 

1 

2 

3 

4 

5 

6 

7 

8 

1 

0.8315 

-0.0029 

0.6053 

-0.0022 

-0.0346 

0.3567 

-0.0347 

0.6192 

2 

0.0299 

0.0059 

-0.0492 

-0.0002 

-0.6041 

-0.1195 

-0.6025 

-0.0040 

3 

0.7498 

-0.0023 

-0.1425 

-0.0007 

0.0614 

-0.9645 

0.0611 

0.2986 

4 

0.0000 

-0.3740 

-0.0012 

0.3563 

-0.0581 

0.0030 

-0.1936 

0.0113 

5 

-0.0014 

-1.0268 

0.0030 

-0.4847 

-0.0434 

-0.0054 

-0.1433 

-0.0304 

6 

0.1430 

0.0067 

0.3071 

0.0047 

0.0000 

-0.3702 

0.0001 

-1.7897 

7 

-0.8114 

0.0024 

0.3292 

0.0031 

0.0010 

-0.5627 

-0.0013 

1.3165 

8 

0.0002 

-1.2687 

-0.0004 

0.9737 

0.3943 

0.0019 

-0.9521 

-0.0072 

9 

1.2653 

-0.0089 

-1.0289 

0.0021 

-0.0009 

0.3151 

0.0164 

0.3747 

10 

-0.0165 

1.0266 

0.0059 

0.0330 

0.6375 

0.0238 

0.6428 

0.0114 


Table 2. Modal Displacements at Control Stations. 


28 







Closed-loop Eigenvalues 

Damping Ratio 

Frequency 

(rad/sec) 

-0.0818 + 0.8139 j 

0.1000 

0.8180 

-0.0036 + 0.8337 j 

0.0043 

0.8338 

-0.0133 + 0.8555 j 

0.0156 

0.8556 

-0.0352 ± 1.1253 j 

0.0313 

1.1258 

-0.2343 ± 1.1112 j 

0.2064 

1.1356 

-0.2059 ± 1.9091 j 

0.1072 

1.9202 

-0.1326 ± 10.7310 j 

0.0124 

10.7319 

-0.2305 ± 14.9345 j 

0.0154 

14.9363 

-0.2526 ± 18.6891 j 

0.0135 

18.6908 

-0.2818 ± 34.0550 j 

0.0083 

34.0561 


Table 3. Closed-loop Eigenvalues With Gain, (! \ . 


Closed-loop Eigenvalues 

Damping Ratio 

Frequency 

(rad/sec) 

-0.0818 ± 0.8139 j 

0.1000 

0.8180 

-0.0830 ± 0.8259 j 

0.1000 

0.8301 

-0.0159 ± 0.8540 j 

0.0186 

0.8541 

-0.1924 ± 0.9840 j 

0.1919 

1.0026 

-0.3145 ± 1.1213 j 

0.2701 

1.1646 

-0.7568 ± 1.9581 j 

0.3605 

2.0993 

-0.2719 ± 10.7216 j 

0.0254 

10.7251 

-1.5302 ± 14.9152 j 

0.1021 

14.9935 

-0.6473 ± 18.7686 j 

0.0345 

18.7798 

-0.4675 ± 34.0439 j 

0.0137 

34.0471 


Table 4. Closed-loop Eigenvalues With Gain, G 2 . 
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Closed-loop Eigenvalues 

Damping Ratio 

Frequency 

(rad/sec) 

-0.0818 + 0.8139 j 

0.1000 

0.8180 

-0.0830 + 0.8259 j 

0.1000 

0.8301 

-0.0857 + 0.8522 j 

0.1000 

0.8565 

-0.3648 + 1.0437 j 

0.3300 

1.1056 

-0.3463 ± 1.1514 j 

0.2880 

1.2024 

-2.5906 + 1.5782 

1.2610 

2.0544 

-2.7128 + 9.5085 j 

0.2744 

9.8879 

-2.0568 + 14.5821 j 

0.1397 

14.7265 

-1.3549 + 18.5311 j 

0.0729 

18.5806 

-0.5327 + 34.0089 j 

0.0157 

34.0131 


Table 5. Closed-loop Eigenvalues With Gain, (i : >. 


Closed-loop Eigenvalues 

Damping Ratio 

Frequency 

(rad/sec) 

-0.0818 + 0.8139 j 

0.1000 

0.8180 

-0.0830 + 0.8259 j 

0.1000 

0.8301 

-0.0857 + 0.8522 j 

0.1000 

0.8565 

-0.3648 + 1.0436 j 

0.3300 

1.1055 

-0.3463 + 1.1514 j 

0.2881 

1.2023 

-2.6015 + 1.5893 

1.2631 

2.0595 

-2.8671 + 9.3545 j 

0.2930 

9.7840 

-3.2694 + 14.3800 j 

0.2217 

14.7470 

-1.8692 + 18.5982 j 

0.1000 

18.6919 

-0.6009 + 34.0359 j 

0.0177 

34.0412 


Table 6 . Closed-loop Eigenvalues With Gain, G' 4 . 
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Closed-loop Eigenvalues 

Damping Ratio 

Frequency 

(rad/sec) 

-0.0818 + 0.8139 j 

0.1000 

0.8180 

-0.0046 ±0.8415 j 

0.0055 

0.8415 

-0.0072 ± 1.1550 j 

0.0062 

1.1550 

-0.0081 ± 1.2455 j 

0.0065 

1.2455 

-0.0997 ± 1.4382 j 

0.0692 

1.4417 

-0.1580 ± 2.0261 j 

0.0778 

2.0322 

-0.0461 ± 10.8179 j 

0.0043 

10.8180 

-0.1222 ± 14.9897 j 

0.0081 

14.9902 

-0.0696 ± 18.7224 j 

0.0037 

18.7225 

-0.0599 ± 34.0745 j 

0.0018 

34.0745 


Table 7. Closed-loop Eigenvalues With Position Gain, G p i and Rate Gain G rl 


Closed-loop Eigenvalues 

Damping Ratio 

Frequency 

(rad/sec) 

-0.0818 ± 0.8139 j 

0.1000 

0.8180 

-0.0830 ± 0.8259 j 

0.1000 

0.8301 

-0.0857 ± 0.8522 j 

0.1000 

0.8565 

-0.1575 ± 1.3598 j 

0.1151 

1.3689 

-0.1682 ± 1.6559 j 

0.1010 

1.6644 

-0.5009 ± 3.7347 j 

0.1329 

3.7682 

-1.0728 ± 10.6741 j 

0.1000 

10.7278 

-1.2597 ± 18.5179 j 

0.0679 

18.5607 

-3.5796 ± 20.9683 j 

0.1683 

21.2717 

-0.7608 ± 34.7506 j 

0.0219 

34.7589 


Table 8. Closed-loop Eigenvalues With Position Gain, G p ± and Rate Gain G , \ 
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Figure 1. Schematic of Phase 2 CEM Structure, With Location of 8 Control Stations. 



